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Abstract 



We present a preliminary study of B 



d Icp)* 

reconstructed in the CP-odd eigenstates K s t:® 



and B 



D® Cp ^K decays, with the D® 



(CP) 

K s u>, in the CP-even eigenstates K + K~, ir + ir~, 
and in the (non-CP) flavor eigenstate K^ir^. Using a sample of about 382 million T(4S) decays 
into BB pairs, collected with the BABAR detector operating at the PEP-II asymmetric-energy B 
Factory at SLAC, we measure the ratios of the branching fractions 



Rcp± = 
and the direct CP asymmetry 



B{B~ -> D° CP± K-) + £(P H 



D° CP± K+) 



[B(B~ -> D°K-)+B{B+ -> D°K+)]/2 



A C p ± 



B(B- 



>D° CP± K- 



B{B~ 



>U CP±^ 



B(B-^D° CP± K~) + B(B+^D° CP± K- 



The results are: 



R CP _ = 0.81 ± O.lO(stat) ± 0.05(syst) 

Pcp+ = 1.07 ± O.lQ(stat) ±0.04(syst) 

A C p- = -0.19±0.12(stat) ±0.02(syst) 

Acp+ = 0.35 ± 0.09(stat) ± 0.05(syst) 
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1 INTRODUCTION 



A measurement of the processes B — > D°K [lj and B — ► D^p^K, where Dq P ^_ indicates the 
CP-even or CP-odd states 1 / \/2(D ± D°) , has been attracting the attention of theorists for the last 
fifteen years [2]. These decay rates are fundamental ingredients in some of the proposed methods 
to extract the 7 angle of the CKM matrix in a theoretically clean way. To this end, one needs to 
measure the two direct CP asymmetries 

_ B(B~ -> D° CP± K-) - B(B+ -> £° P± if+) 
CP± -► Dg^ir-) + £(P+ D° CP± K+) ' 1 j 

and the two ratios of charge-averaged branching fractions in D° decays to CP eigenstates 

g(g- -» D° C p±K-) + B(g+ -^g P± iT+) 
CP± [B(B~ D°/jf-) + £(P+ -» /2 ' U 

In fact, 7 is constrained by the following set of equations in the three unknowns 7, 8, r: 

Rcp± = 1 + r 2 ±2r cos 8 cos 7 (3) 
±2rsin<5sin7 

Acp± = 1^210 r — ' ( 4 ) 

1 + r z ± 2r cos cos 7 

where r = \A(B~ -> B°K~)/A(B- -> 0(0.1) is the magnitude of the ratio of the 

amplitudes for B~ — > D°K~ and P~ — > D°K~ and <5 the difference between their strong phases. 
The asymmetries ^4cp±, in addition to being ingredients for the extraction of 7, are of special 
relevance because they would indicate, if significantly different from zero, direct CP violation in 
charged B decays. To measure Rcp± and Aqp± we reconstruct B — ► Dq P± K and P» — ► Z)gp_|_7r 
decays with the -D°p_|_ decaying to two CP-odd and two CP-even eigenstates, and B — ► D°K~ 
and P~ — > D°ir~ decays with Z)° decaying to one non-CP state. Previous measurements of these 
quantities were performed by BABAR [3 J and Belle [1]. We update the result by BABL4R from 211 fb _1 
to 348 fb -1 of data. Compared to the previous analysis, the current study does not include the 
decay mode D° — ► Kg(f>, since it is going to be explored by the Dalitz analysis method using 
B~ ^DK~ , D^KgK + K~ decays. Dropping D° — ► K®<j> allows to combine the results of both 
measurements in the future. We also express the CP-sensitive observables in terms of three Dalitz 
related independent quantities: 

R C p+{lT A CP+ ) -Rcp-{1^A CP ~) 
x± = 4 , (5) 

r* = X l + yl= RcP+ + * CP --\ (6) 

where the Cartesian coordinates x± = r cos(5±7) and y± = r sin(5±7) are the same CP parameters 
as were measured by the BABAR Collaboration using B~^DK~ , D^Kgir~-K + decays [5]. We 
reduce the systematic uncertainties from D° branching fractions and reconstruction efficiencies of 
different D° modes by measuring the double branching fraction ratios 

R KI " 
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rather than the quantities Rcp±- Here, 



Kj, = B{B~ -» D° CP± K-) + B(B+ -» D Q cp ± K+) 
,( ' ,,± - B(B--23g p± 7r-)+B(B+ ^^ P± 7T+) ' 



n CP±— . r>0 _-\ , nfDA- . n(l _4-\ ' \°) 

and 



R K/« = g(i3-^ J P ^-)+g(i3 + ^g ^ + ) 

-» DOvr-) + £(P+ -► ' 1 ^ 

i?± and Rcp± are equivalent discarding a term of the order of ~ 0.01, which will be accounted for 
by assigning a systematic uncertainty when quoting the result in terms of Rcp±- 



2 THE BABAR DETECTOR AND DATASET 

This measurement uses 348 fb^ 1 of data taken at the T(4S) resonance by the BABAR detector 
with the PEP-II asymmetric B factory. The BABAR detector is described in detail elsewhere [6]. 
Tracking of charged particles is provided by a five-layer silicon vertex tracker (SVT) and a 40- 
layer drift chamber (DCH). A ring-imaging Cherenkov detector (DIRC) provides improved particle 
identification (PID). An electromagnetic calorimeter (EMC), comprised of Csl crystals, is used to 
identify electrons and photons. These systems are mounted inside a 1.5 T solenoidal magnetic field. 
The instrumented flux return of the magnet allows discrimination of muons from other particles. 
We use a GEANT4-based Monte Carlo (MC) simulation [7] to model the response of the detector, 
taking into account the varying accelerator and detector conditions. 



3 EVENT SELECTION 

We reconstruct B — > D°h~ decays where the prompt track h~ is a kaon or a pion. Candidates 
for D° are reconstructed in the CP-even eigenstates ir~ir + and K~K + , in the CP-odd eigenstates 
K°tt°, K° s uj and in the non-CP flavor eigenstate K tt + . and u candidates are selected in the 
7r + 7r~ and 7r + 7r~7r° channels, respectively. 

PID information from the DCH and, when available, from the DIRC must be consistent with 
the kaon hypothesis for the K meson candidate in all D° modes and with the pion hypothesis for 
the tt^ meson candidates in the D°^ir~ir + mode. For the prompt track to be identified as a pion 
or a kaon, we require at least five Cherenkov photons to be detected to ensure a good measurement 
of the Cherenkov angle. We reject a candidate track if its Cherenkov angle is not within 4er of the 
expected value for either a kaon or pion mass hypothesis. We also reject candidate tracks that are 
identified as electrons by the DCH and the EMC or as muons by the DCH and the muon system. 

Photon candidates are clusters in the EMC that are not matched to any charged track, have 
a raw energy greater than 30MeV and lateral shower shape consistent with the expected pattern 
of energy deposit from an electromagnetic shower. Photon pairs with invariant mass within the 
range 115-150 MeV/c 2 (~3cr) and total energy greater than 200 MeV are considered 7r° candidates. 
To improve the momentum resolution, the 7r° candidates are kinematically fit with their mass 
constrained to the nominal tt° mass [8]. 

Neutral kaons are reconstructed from pairs of oppositely charged tracks with invariant mass 
within 10 MeV/ c 2 from the nominal K® mass [8]. We also require the ratio between the flight 

7 The double branching fraction ratios, in the approximation A(B + — > Dq P ±tt + ) ~ A(B~ — > Dcp±7r~) ~ 
-i A(B~ — > D°ir~), are equivalent to Rcp±, discarding a term rB\V ua V c d/V u dVcs\ ~ 0.01. 
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distance in the plane transverse to the beam direction and its expected uncertainty to be greater 
than 2. 

The invariant mass of a D° candidate must agree within 2.5a of its mass resolution to the 
nominal D° mass [8j. The D° mass resolution is about 7.5 MeV/c 2 in the Kir, K + K~ and 7r + 7r~ 
modes, and about 21 MeV/c 2 and 9 MeV/c 2 in the K^ir and K®uj modes, respectively. Selected 
D° candidates are fit with a constraint to the nominal D° mass. 

We reconstruct B meson candidates by combining a D° candidate with a charged track h~ . For 
the K~ir + mode, the charge of the track h~ must match the one of the kaon from the D° decay. 
We select B meson candidates by using two kinematically independent variables: the beam-energy- 
substituted mass 

m ES = ^(Ef/2 + Pl . PB r/E?-pl 

and the energy difference 

AE = E* B -E*/2, 

where the subscripts i and B refer to the initial e + e~ system and the B candidate, respectively, and 
the asterisk denotes the beams center-of-mass (CM) frame. The mES distributions for B — ► D°h~ 
signal events are Gaussian distributions centered at the B mass with a width of 2.6MeV/c 2 , which 
does not depend on the decay mode or on the nature of the prompt track. In contrast, the AE 
distributions depend on the mass assigned to the prompt track. We evaluate AE with the kaon 
mass hypothesis so that the distributions are centered near zero for B~ — ► D°K~ events and shifted 
on average by approximately 50MeV/c to the positive direction for B — > D°-k~ events. The AE 
resolution depends on the momentum resolutions of the D° meson and the prompt track h~~, and 
is typically 16 MeV for all the D° decay modes. All B candidates are required to have uies within 
2.5a of the mean value and AE in the range —0.15 < AE < 0.20 GeV. 

To reduce background from continuum production of light quarks, we construct a Fisher discrim- 
inant based on the following four quantities: (i) The Legendre polynomials, a set of momentum- 
weighted sums of the tracks and neutrals not associated with the reconstructed candidate, i.e. 
coming from the rest of the event (ROE): 

ROE 

Lj = x l cos (W, ( 10 ) 

i 

where 9* is the CM angle between p* and the thrust axis T B of the B candidate. We have 
considered only the Lq, L2 pair, since previous studies have shown that adding other Lj to the set 
of discriminating variables does not improve the sensitivity. In particular we use the ratio L2/L0; 
(ii) i? ROE , the ratio of the Fox- Wolfram moments Hf'^/H^ ^', computed using tracks and photons 
in the ROE. i^ ROE is defined as [9]: 

^r^E^P^cos^), (ii) 

i,j vis 

where P; is a Legendre polynomial, 9ij is the opening angle between p* and p*, and E* is is the 
total visible energy of the event, (iii) |cos(p^,z)| is the cosine of the angle of the B candidate 
momentum with respect to the beam (z) axis, (iv) | cos(T B , z)\ is the cosine of the angle of the 
B candidate thrust axis with respect to the z axis. A cut on the value of the Fisher discriminant 
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rejects more than 75% of the continuum background while retaining about 85% of the signal in all 
modes. 

Another source of background is related to BB events. Its main contributions come from the 
processes B^D*h (h = tt,K) and B~^D°p~ mis-reconstructed as B — > D°/i~candidates. For 
D° — > K~K + , D° — > 7r _ 7r + , D° — > K^tt and D° — » K®uj decays, there are peaking backgrounds 
caused by B mesons decaying into the same final state particles. The peaking backgrounds have 
AE and tties distributions similar to the signal. Their yields are estimated from the D° invariant 
mass sideband data and are taken into acount in the fit. 

When reconstructing B candidates, it is possible that more than one combination satisfies 
the selection criteria in the same event. In order to select only one candidate per event, we 
define a criterion that allows to identify the combination with the largest probability to be a true 
signal B~ —> D°h~ decay. The D° invariant mass and the energy-substituted mass are chosen as 
discriminating quantities in all the channels. When D° decays into the CP-odd channels we also 
include the invariant masses of the uj and ir° candidates. These variables are combined in a x 2 
function whose minimization defines the best candidate choice. In the end, the fraction of rejected 
background candidates in the selected samples is 2% in the Ktt, less than 1% for KK and tt + tt~ 
modes, while it is about 5% in the K^n mode and 8% in the K®u). 

The total reconstruction efficiencies, based on simulated signal events, are about 35%(K~ir + ), 
32%(K-K+), 33%(vr-7r+), 20%(i^°) and 8%(K° s uj). 



4 FIT PROCEDURE 

We determine the signal and background yields for each D° decay mode from a two-dimensional 
extended unbinned maximum-likelihood fit to the selected data events determines the signal and 
background yields. The input variables to the fit are AE and a particle identification probability 
for the prompt track based on the Cherenkov angle 9c, the momentum p and the polar angle 9 of 
the track. The extended likelihood function C for the selected sample is given by the product of 
the probabilities for each individual candidate and a Poisson factor: 

p -N' (N /\N 

i=l 

The probability V% for a candidate in the event i is the sum of the signal and background terms: 

Vi{AE,0 ) = ^vP^ + ^Vr K + (13) 
N qqM p&fr) + N m(K) v qq{K) + 



N' 1 N 
N BB(-k) v bb{k) N BB(K) v bb{K) 

JSfi ' N' 1 

N X 1 X 2 K pX!X 2 K 



N' 

where N' = N D o 7T +N D o K +N q g^ ) +N q ^ K) +N BS{7T) +N BS ^ K) +N Xl x 2 K- Each addendum on the 
right-hand side of equation (|13p is the product of two different terms. The ratio Nj/N' (J = 
D°ir, D°K, ...) represents the probability to choose a candidate of type J after the selection criteria 
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are applied; the term V( is the probability to measure the particular set of physical quantities 
{AE, Oc}i in the i th event, once the candidate of type J has been selected: 

vi = vi Eti v J 6c>i . (14) 

The AE distribution for B~ — ► D°K~ signal is parameterized with a double Gaussian function. 
The mean and width of the narrow Gaussian are denoted in the following with fi(D°K) and <j(D°K). 
The B~ — » D°7r~ AE probability density function (PDF) would be the same as the B — » D°K~ 
one, if the prompt track had been assigned the pion mass. Since AE is computed by assigning the 
kaon mass, it is shifted by a quantity 

AE sUit ( 7 , \p\) = 7 (yJm 2 K + \0\2 - vX + IPl 2 ) 

which depends on the momentum p of the prompt track in the lab frame. 7 is the Lorentz parameter 
characterizing the boost of the center of mass frame relative to the lab frame. Therefore we 
parameterize the B~ — ► D°tt~ AE shape with a double Gaussian whose mean is computed, event- 
by-event, as ij>(D it) = fi(D°K) + Ai£ s hift(7, |p|)> anci whose width is the same as for the B~ — > 
D°K~ signal component. The fraction of the wide component of the signal shape, its offset from the 
narrow component and the ratio between its width and the width of the narrow component are fixed 
using the mode-dependent numbers obtained from the MC simulation. The AE distributions for the 
continuum background are parameterized with a first order polynomial. The AE distribution for 
the BB background is empirically parametrized with a "Crystal-Ball" lineshape [10J: a Gaussian 
with an exponential tail at higher AE values. The parameters of the background shapes are 
determined from MC simulated events and are fixed in the fit. 

The particle identification PDF is obtained from MC simulation. Its parametrization is per- 
formed by means of a double Gaussian distribution as a function of 9q^~, which is the difference 
between the measured Cherenkov angle and its expected value for a given mass hypothesis, divided 
by the expected error. 

We independently fit five samples corresponding to each of the five D° decay modes under study. 
The fit simultaneously evaluates separate likelihood functions for B + and B~ categories. In the fit 
the free parameters are D K and D°7r signal yield asymmetries, total number of signal events in 
D°tt (A^oJ, ratio R K/n = N D o K /N D07n eight background yields: N q ^ K) , N bS(7t) , N bE{k) 

(one for each charge, i.e. 4 x 2 = 8), and two parameters of the AE signal shape (shared between 
positive and negative samples). The number of peaking background events Nx 1 x 2 K is fixed to the 
values obtained from the study using D° mass sidebands. We assume no charge asymmetry in the 
peaking background (a small systematic error due to this assumption is considered later). 

5 PHYSICS RESULTS AND SYSTEMATIC STUDIES 

The results of the fit are summarized in Table [TJ Figure Q] shows the distributions of AE and 9c for 
the K~~tt + , CP-even and CP-odd modes. The projections of the likelihood fits are overlaid on the 
plots. On Figure [2] we show the AE projections produced with a kaon identification requirement 
applied to the prompt track. Hence the B~ — ► D°K~ signals become prominently visible on the 
plots, while B — > D°ir~ contributions significantly decrease. 

The double ratios Rcp± are computed by calculating a weighted mean of the ratios R^u for 
CP-even and CP-odd modes and dividing it by Rk/tt f° r the non-CP mode. Correction factors 
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Table 1: B — > and P — > D°7r signal event yields obtained from the fit to the data. All 

values are preliminary. 



D {) mode 


N(B -» P> u vr) 


N(B D u i^T) 


N(B+ -> D°K+) 


N(B~ -» £ U K~) 


K-7T+ 


24965 ± 169 


1859 ± 52 


951 ± 36 


909 ± 35 


K-K+ 


2412 ± 54 


189 ± 20 


61 ±11 


128 ±16 


7T~7T + 


876 ± 38 


73 ± 15 


24 ±9 


49 ±12 


K u s tt u 


2967 ± 69 


184 ± 25 


107 ±19 


77 ± 16 




1023 ± 44 


59 ± 14 


36 ± 12 


23 ±9 



(ranging from 1.006 to 1.027 depending on the D° mode) that account for small differences in 
the efficiency between the B~ — > D°K~ and B~ — > D°ir~ selections are taken into account. An 
additional factor is applied to the results in the D° — > K^uj mode to correct for a dilution due 
to the S-wave non-resonant contribution. These corrections were estimated using a fit to the u 

helicity angle in the selected data events and found to be 1.10 ± 0.11 for A C p_ and 0.98 ± 0.02 for 

Rqp^_. The uncertainties in the correction factors are included in the systematic errors (±0.006 
and ±0.008 for Rcp- and Acp-, respectively). The results for each mode separately and combined 
by CP-even and CP-odd categories are listed in Table [2j 

Table 2: Measured double branching fraction ratios Rcp± and CP asymmetries Aqp± for different 
D° decay modes. In the combined results, the first error is statistical, the second is systematic. 
For individual modes, only statistical errors are shown. All values are preliminary. 



Z) u decay mode 


Rcp 


A C p 


K~K+ 


1.05 ±0.11 


0.36 ±0.10 




1.13 ±0.22 


0.33 ±0.20 


CP-even combined 


1.07 ±0.10 ±0.04 


0.35 ±0.09 ±0.05 




0.84 ±0.11 


-0.16 ±0.13 


K° s lo 


0.75 ±0.18 


-0.24 ±0.26 


CP-odd combined 


0.81 ±0.10 ±0.05 


-0.19 ±0.12 ±0.02 



Systematic uncertainties in the double ratios Rcp± and in the CP asymmetries Aqp± arise pri- 
marily from uncertainties in signal yields due to the estimate of the peaking backgrounds (±0.03 for 
Rcp+ and ±0.05 for Rcp~) and from the imperfect knowledge of the PDF shapes. The parameters 
of the PDFs that are fixed in the nominal fit are varied by ±ler and the observed difference in the 
parameters Rk/-ki signal yield asymmetries and N D o K is taken as a systematic uncertainty (±0.003 
for A C p+, ±0.002 for A CP -, ±0.010 for R C p+ and ±0.007 for R C p-)- Possible CP asymmetries 
up to 20% in the peaking backgrounds are also taken into account (±0.04 for Aqp+)- An estimate 
of the intrinsic detector charge bias due to acceptance, tracking, and particle identification effi- 
ciency has been obtained from the weighted average of the measured asymmetries in the processes 
B~ — >D°[— >K~ir + ]h~ and B~—*-D^p,ir~, where CP violation is expected to be negligible. This 
asymmetry estimate (±0.02) has been added in quadrature to the total systematic uncertainties on 
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the CP asymmetries Acp± (this is a correlated part of the systematics for Acp+ and Acp- ) • The 
accuracy in the equivalence between R± and Rcp± is evaluated to be ±0.03 for Rcp+ and ±0.02 
for Rcp- (these uncertainties are correlated). 

6 SUMMARY 

In conclusion, we reconstruct B~ — > D°K~ decays with D° mesons decaying to non-CP K^ir , 
CP-even K + K~ and 7r + 7r~ and CP-odd K®7r° and K®u> eigenstates. We have measured the CP 
asymmetries A C p+ = 0.35 ± 0.09(stat) ± 0.05(syst) and A CP - = -0.19 ± 0.12(stat) ± 0.02(syst). 
Our result for Acp+ is 3Aa away from 0. This constitutes the first evidence for direct CP violation 
in B~ — ► D°K~ decays. The double ratios of branching fractions are measured to be Rcp+ = 
1.07 ± O.lO(stat) ± 0.04(syst) and R C p- = 0.81 ± O.lO(stat) ± 0.05(syst). 

The corresponding values of x± and r 2 B are extracted using equations 5 and 6, separately 
propagating correlated and uncorrelated errors on Acp± and Rcp±- We obtain x + = —0.065 ± 
0.047(stat) ± 0.020(syst), x_ = 0.199 ± 0.052(stat) ± 0.020(syst), r| = -0.060 ± 0.070(stat) ± 
0.039(syst). 

The results obtained in this analysis are statistically in agreement with the previous measure- 
ments as demonstrated in Table [3l All results presented in this document are preliminary. 

Table 3: Comparison of the preliminary results of this analysis to the previous measurements by 
BABAR [3] and Belle [I]. The decay mode D° — > Kg(f>, used in the previous analyses, is not included 
in the present measurement. 



Parameter 


Present analysis 


BABAR (2006) [3j 


Belle (2006) [4J 


Rcp- 


0.81 ±0.10 ±0.05 


0.86 ±0.10 ±0.05 


1.17 ±0.14 ±0.14 


Rcp+ 


1.07 ±0.10 ±0.04 


0.90 ±0.12 ±0.04 


1.13 ±0.16 ±0.08 


A C p- 


-0.19 ±0.12 ±0.02 


-0.06 ±0.13 ±0.04 


-0.12 ±0.14 ±0.05 


Acp+ 


0.35 ±0.09 ±0.05 


0.35 ±0.13 ±0.04 


0.06 ±0.14 ±0.05 
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Figure 1: AE (left) and 9c (right) distributions of selected B~ — * D°h~ events. The blue line 
represents the projection of the likelihood in the plotted variable. The red line represents the 
B~ — > D°K~ component. In the left hand plots, the green and light blue lines indicate BB and 
continuum backgrounds, respectively. The brown line refers to the peaking background (when 
present). 
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Figure 2: AE distributions of B~ — > D°K~ signal enhanced B~ — > D°h~ events. The blue line 
represents the projection of the likelihood, the red line indicates the B~ — > D®K~ component, the 
green line shows the total background contribution. The remaining B~ — ► D°tt~ signal is visible 
as a small shoulder on the right hand side of the B~ — > D°K~ signal peak. 
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